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Semiconductor nanowhiskers made of III-V compounds exhibit great potential for technological 
applications. Controlling the growth conditions, such as temperature and diameter, it is possible to 
alternate between zinc blend and wurtzite crystalline phases, giving origin to the so called polytyp- 
ism. This effect has great influence in the electronic and optical properties of the system, generating 
new forms of confinement to the carriers. A theoretical model capable to accurately describe elec- 
tronic and optical properties in these polytypical nanostructures can be used to study and develop 
new kinds of nanodevices. In this study, we present the development of a wurtzite/zincblend poly- 
typical model to calculate the electronic band structure of nanowhiskers based on group theory 
concepts and the k-p method. Although the interest is in polytypical superlattices, the proposed 
model was applied to a single quantum well of InP to extract the physics of the wurtzite/zincblend 
polytypism. By the analysis of our results, some trends can be predicted: spatial carriers' sep- 
aration, predominance of perpendicular polarization (xy plane) in the luminescence spectra and 
interband transition blueshifts with strain. A possible range of values for the WZ InP spontaneous 
polarization was suggested. 



I. INTRODUCTION 

Low dimensional semiconductor structures exhibit dif- 
ferent characteristics, ruled by the size and morphology 
of the system. Recently, there is an increasing interest 
in nanowhiskers (NWs), also known as nanowires. These 
are one dimensional nanostructures grown perpendicular 
to the surface of the substrate, usually by the vapor- 
liquid-crystal (VLC) method. The technological applica- 
tions of NWs, including biological and chemical nanosen- 
sors [l|-|3], lasers [4|, light emission diodes [5] and field 
effect transistors [6], can be used in a large variety of 
fields. 

The first register in the literature of whiskers was made 
by Wagner and Ellis 0] in 1964. In this classic study, it 
is demonstrated the vertical growth of a Si whisker in 
the [111] direction, activated by droplets of Au, using 
the VLC method. The radius of the structure is approx- 
imately the same as the catalist droplet of Au and the 
vertical size depends on the ratio and time of the com- 
pound's deposition on the substrate. Although the VLC 
method is the most common, other methods like vapor- 
phase epitaxy (VPE), molecular-beam epitaxy (MBE) 
and magnetron deposition (MD) are also applied for the 
NWs growth. 

In III-V compound NWs (e. g., arsenides and phos- 
phides), a surprising characteristic is the predominance 
of the wurtzite (WZ) phase in the structure. Exception 
made for the nitrides, the stable crystalline structure of 
III-V compounds in the bulk form is zincblend (ZB). Al- 
though the difference between the small formation energy 
of the two phases are small, approximately of 20 meV 
per pair of atoms at zero pressure, high pressures would 
be necessary to obtain the WZ phase in the bulk form. 
However, reducing the dimensions of the system to the 



nanoscale level, such as in these NWs, the WZ phase be- 
comes more stable. This stability is due to the smaller 
surface energy of lateral faces compared to the cubic crys- 
tal. An extensive summary of NWs growth, properties 
and applications was made by Dubrovskii et at. [8|. 

Controlling the growth conditions, such as tempera- 
ture and diameter of the NW, it is possible to create 
different regions composed of ZB and WZ structures 
16]. The mixture of both crystalline phases in the same 
nanostructure is called polytypism. Such characteristic 
directly affects the electronic and optical properties of 
NWs. The detailed study of polytypism in III-V semi- 
conductor NWs is fundamental to the development of 
novel functional nanodevices with a variety of features. 

The theoretical tool to calculate the electronic band 
structure of polytypical NWs used in this paper is the 
k-p method. Although the formulation of this method 
has already been done for ZB and WZ crystal structures 
in the bulk form I17H21] and in superlattices and het- 
erostructures [22.- 26], it was never applied to a polytyp- 
ical case. 

Deeply studying the core of the k-p formulation for 
both crystal structures and relying on the symmetry con- 
nection of the polytypical interface presented in the paper 
of Murayama e Nakayama [i^ , it was possible to describe 
the top valence bands and lower conduction bands of ZB 
and WZ in the same Hamiltonian matrix. The envelope 
function scheme was then applied to obtain the variation 
of the parameters along the growth direction, describing 
the different regions in a NW, thus completing the poly- 
typical model. The effects of strain, spontaneous polar- 
ization and piezoelectric polarization are also included in 
the model. 

In order to test the model, it was applied to a polytyp- 
ical WZ/ZB/WZ quantum well of InP. Although a real 
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NW is composed by several regions of WZ and ZB, the 
physical trends of the polytypical interface can be ex- 
tracted from a single quantum well system. We choose 
the InP compound basically for two reasons: the small 
spin-orbit energy makes easier to fit the matrix parame- 
ters to the effective mass values given in the paper of De 
and Pryor [28] for the WZ polytype and also because the 
InP NWs can be found in a great number of studies in 
the present literature [ol, l29l435| . 

The present paper is divided as follows: In section II we 
discuss the symmetry of ZB and WZ crystal structures 
and analyze how the irreducible representations of the 
energy bands are connected in the polytypical interface. 
Section III describes the Hamiltonian terms for the poly- 
typical model. The results, and their discussion, of InP 
WZ/ZB/WZ single well are found in section IV. Finally, 
in section V, we draw our conclusions. 





ZB 



WZ 



Figure 1. ZB (left) and WZ (right) structures and their stack- 
ing sequence. The ZB is presented in the [111] direction. In 
this direction, the unit cell is twice as large as the usual ZB 
unit cell |27|]. 



Irreducible representations at the polytypical 
interface 



II. SYMMETRY ANALYSIS 



A. Zincblend and wurtzite structures 



Our formulation relies on group theory concepts and 
therefore it is necessary to understand the symmetry of 
the two crystal structures considered in the polytypical 
NWs. The ZB structure belongs to the class of symmor- 
phic space groups and has the symmetry as its factor 
group. The number of atoms in the primitive unit cell is 
two. Unlike ZB, the WZ structure belongs to the class 
of nonsymmorphic space groups and its factor group is 
isomorphic to Csy. The classes of symmetry operations 
C2, Cq and (jy are followed by a translation of c/2 in the 
[0001] direction. WZ has four atoms in the primitive unit 
cell. Comparing the factor groups one can notice that 
Cqv is less symmetric than T^. In the k-p framework, this 
lower symmetry decreases the number of irreducible rep- 
resentations (IRs) in the group consequently increasing 
the number of interactions in the Hamiltonian. A good 
description of the concepts of space group symmetry can 
be found in Ref. [36]. 

In polytypical NWs, the common growth direction is 
the ZB [111], which exhibit a noticeable similarity to WZ 
[0001]. Actually, analyzing both crystal structures in 
these directions, one can describe then as stacked hexag- 
onal layers. The ZB has three layers in the stacking se- 
quence (ABCABC) while WZ has only two (ABAB) as 
shown in Figure [TJ The crystal structure alternation oc- 
curs when a stacking fault happens in WZ, leading to a 
single ZB segment, or when two twin planes appear in 
ZB, originating a single WZ segment |ll| . 



An important issue of the model is how to connect the 
energy levels at the polytypical interface depending on 
the symmetry they assume. Based on the scheme pre- 
sented by Murayama and Nakayama [27] of single group 
IRs at the WZ/ZB interface the symmetry of the energy 
bands can be chose. The same scheme was constructed 
by De and Pryor [isl] for the double group IRs with the 
inclusion of the spin-orbit coupling. 



A [0001] 




Figure 2. Usual a) ZB and (b) WZ first Brillouin zones and 
their respective high symmetry points. The arrows (red line) 
represent the growth directions in NWs. For ZB, the [111] 
direction is directed towards the L-point. 

Since WZ has twice more atoms in the primitive unit 
cell than ZB, the number of energy bands in the F-point 
is also twice as large. Considering the sp^ hybridization, 
without spin, ZB has 8 energy bands while WZ has 16. 
However, in the [111] direction, the ZB unit cell is two 
times larger than the usual face-centered cubic (FCC) 
unit cell [27|. In the IRs scheme mentioned above, the 
presence of energy bands with L symmetry takes into 
account the mismatch in the number of atoms for the 
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usual unit cells. The reason for the appearance of the L 
symmetry is the fact that ZB [111] growth direction is 
directed towards the L-point, as displayed in Figure [2^, 
hence this point is mapped out in the F-point. Figure 
[SJd, displays the first Brillouin zone (FEZ) for the WZ 
structure. 
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Figure 3. The subset of IRs considered in this formulation 
with and without the spin-orbit (SO) coupling. The numbers 
in parentheses are the degeneracy of the IRs. The notation 
for the IRs follow Refs. [23, H 



Among all the IRs presented in Refs. [27|, |28| we con- 
sidered only a small subset. Displayed in Figure [H this 
subset comprises the lower conduction band and the top 
three valence bands, which belong to the F-point in both 
structures. The price that is paid in considering only 
a small subset is the accuracy of the Hamiltonian for a 
fraction of the FEZ, approximately 10-20%. The basis 
states for the IRs of the considered bands are presented 
in equations ([T]) and (j2j) for ZE and WZ, respectively. 
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Although the IRs belong to different symmetry groups 
{Td and Ce^), the basis states transform as the usual 
x^y^z cartesian coordinates for the valence bands and the 
scalar x^+^^ + z^ for the conduction band in both crystal 
structures. This information is crucial to represent WZ 
and ZE with the same Hamiltonian and is the essential 
insight of our formulation. 



III. THEORETICAL MODEL 



A. k p Hamiltonian 



In order to develop our k-p Hamiltonian [375 conve- 
nient to describe the ZE structure in a coordinate system 
that has the z axis parallel to the growth direction. This 
coordinate system is the primed one presented in Figure 
m Even though the choice of the coordinate system is 
arbitrary, it alters the k-p Hamiltonian. For example, in 
the unprimed coordinate system the kz direction is di- 
rected towards the X-point but in the primed coordinate 
system it is directed towards the L-point. Thus our ex- 
pectation is that an anisotropy in the ZE Hamiltonian 
between the kx and kz directions will occur since they 
are not reaching equivalent points in the the reciprocal 
space anymore. 
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Figure 4. (a) ZB conventional unit cell with two different 
coordinate systems, (b) WZ conventional unit cell with its 
common coordinate system. The [111] growth direction for 
ZB structure passes along the main diagonal of the cube and 
is represented in the primed coordinate system. 



Considering the single group formulation of the k-p 
method, the choice of the coordinate system defines the 
symmetry operation matrices used to derive the momen- 
tum matrix elements. It is necessary to recalculate the 
Hamiltonian terms for the ZE structure. However, the 
energy bands we consider here are exactly the ones cus- 
tomary used in the ZE [001] k-p Hamiltonian. Instead 
of recalculating the terms for ZE [111] k-p Hamiltonian 
it is possible, and also useful, to apply a basis rotation 
to the ZE [001] matrix. This rotation procedure is well 
described in the paper of Park and Chuang [38] . 

The basis set for both crystal structures in the primed 
coordinate system (the prime will be dropped out of the 
notation from now on and will be used only when it is 
necessary) is given by 
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In a first approximation, tiie interband interaction is 
not taken into account explicitly here, thus the conduc- 
tion band is a single band model for spin-up and spin- 
down reading as 
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(4) 



where Eg is the band gap, £^0 is the energy reference 
at ^ = and mH, are the electron effective masses 
parallel and perpendicular do the z axis, respectively. For 
the ZB structure, however, the electron effective masses 
are equal. 

The Hamiltonian for WZ and ZB valence band is given 

by 
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and the matrix terms are defined as 
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Aikl^A2 {kl - 
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V2A3 



A^kt 



(6) 



coordinate system, Ai{i = l,...,6,z) are the holes ef- 
fective mass parameters, Ai is the crystal field splitting 
energy in WZ, A2,3 are the spin-orbit coupling splitting 
energies and k± = kx ^ iky. 

It is important to notice that the parameter A^ appears 
in the matrix elements to regain the original isotropic 
symmetry of the ZB band structure in the new coordi- 
nate system. In the regions of WZ crystal structure, this 
parameter is zero and the matrix is exactly the canonical 
in use for WZ crystals. 

Although this is not the usual way to describe ZB crys- 
tals, all the parameters in the matrix can be related to 
the familiar = 1, 2, 3) and A so as shown above 



Ai = 
A2 = A3 = 

Ai -71 - 473 
A2 -71 + 273 
A3 673 

-373 
^5 -72 - 273 
A6 = -V2(272 H 
Az =^2- 73 



73) 



(7) 



where ka{a = x, z) are the wave vectors in the primed 



One may also argue that this formulation is very sim- 
ilar to the WZ phase. The insight here is to consider ZB 
as a WZ structure without the crystal field splitting en- 
ergy. Since the WZ structure is less symmetric than ZB, 
as mentioned in section II. A, it is possible to represent 
the ZB parameters with the WZ ones. 

The resulting valence band structures for bulk WZ and 
ZB InP using matrix ^ are shown in Figure [H The pres- 
ence of the crystal field in the WZ structure creates three 
distinct two-fold degenerate bands whereas in ZB there 
is a four-fold and a two-fold degenerate set of bands. Ad- 
ditionally, the anisotropy between kx and kx is evident 
in both crystal structures. In WZ it is due to the differ- 
ent symmetry properties of the xy-plane and the z axis 
but in ZB it is because kx and kz directions do not reach 
equivalent points in the reciprocal space. Since the con- 
duction band is a parabolic model, we did not present its 
dispersion relation. The ZB parameters were obtained 
from Ref. [39] and the WZ parameters were derived us- 
ing the effective masses presented in Ref. [28]. These 
parameters can be found in Table |T1 

Following Chuang and Chang [40] notation, the valence 
energy bands for WZ are named after the composition of 
the states at /c = 0. HH (heavy hole) states are composed 
only by |ci) or |c4), LH (light hole) states are composed 
mainly of |c2) or IC5) and CH (crystal- field split-off hole) 
are composed mainly of IC3) or Ice). For the ZB structure 
the common identification of the valence energy bands 
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Figure 5. Valence band structure for bulk (a) WZ and (b) ZB 
in the primed coordinate system. The usual identification of 
the bands was used for ZB while in WZ it was necessary to 
analyze the composition of the states in the F-point. We can 
see the anisotropy between kz and kx in WZ and also in ZB 
because in the new coordinate system the x and z axes do 
not reach equivalent points in the reciprocal space. The top 
of the valence band in both crystal structures was chosen to 
be at zero. 



was used. The four-fold degenerate bands at /c = are 
HH and LH and the lower two-fold degenerate band is 
SO (split-off hole). 



Table L InP parameters used in the calculations. 



Parameter 



ZB InP WZ InP 



Lattice constant (A) 

a 5.8697 4.1505 

c - 6.7777 
Energy parameters (eV) 

Eg 1.4236 1.474 

Ai 0.303 

As = As 0.036 0.036 
Conduction band effective masses 

ml/mo 0.0795 0.105 

m^/mo 0.0795 0.088 
Valence band effective mass 
parameters (units of 

Ai -13.4800 -10.7156 

A2 -0.8800 -0.8299 

As 12.6000 9.9301 

A4 -6.3000 -5.2933 

-5.8000 5.0000 

Aq -7.4953 1.5000 

A^ -0.5000 



where a^w and ac± are the conduction band deformation 
potentials parallel and perpendicular to the z axis, re- 
spectively. In the ZB structure they have the same value. 
The valence band strain Hamiltonian is 
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(9) 



and the matrix terms are 



B. Strain 

The strain Hamiltonian can be obtained using the same 
basis rotation applied to the k-p matrix [38|. Similarly, 



the conduction and valence band are decoupled. 

For the conduction band, the strain effect is given by 



Ece = CLc\\£zz + CLcl. {^xx + ^yy) 



(8) 



F, = {Di + Ds) 6,, + iD2 + D^) {s. 



-yy 



-^5 ^ {,^xx 



£yy) + ^2iexy 



-yyJ 



(10) 



where I^^'s are the valence band deformation potentials 
and Eij is the strain tensor. 

In the same way as the A^ parameter appears in k-p 
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matrix, some extra deformation potential terms were ap- 
pears to use the same strain Hamiltonian for both crys- 
tal structures. The deformation potential Dr^ was split 
in two parts because the strain tensor e^y is not present 

in the ZB structure. For WZ, D^^^ = Df\ Also, the 
deformation potential takes into account the non existing 
term £xx — £yy in the WZ structure. 

The deformation potentials I^^'s are related to the ZB 



where ao is the lattice constant of the substrate. 

In the unprimed coordinate system, the strain tensor 
elements assume the form 
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Considering biaxial strain, the elements of the strain 
tensor can be obtained in both coordinate systems for 
ZB and WZ. Although it is convenient to describe the 
Hamiltonian terms in the primed coordinate system, it 
is also useful to describe the strain tensor elements in 
the unprimed coordinate system. They will be used to 
construct the piezoelectric polarization in section III.C. 
The prime will be reintegrated in the notation to avoid 
confusion and the upper scripts z and w denotes the ZB 
and WZ structures, respectively. 

For the primed coordinate system, the elements of the 
strain tensor are given by 



-xy 



1 / _j_ 

"3 r + ct(iii) 



(17) 



The quantity a^^^^^ is given by 
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Comparing the expression (fT3|) with (p!4|) it is possible 

(z) (z) 

to obtain effective values for C^^ and C33 for ZB in the 
primed coordinate system. The effective values are: 
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Thus, we have a single set of expressions to describe 
biaxial strain in the primed coordinate system for ZB and 
WZ crystal structures: 



ao — <2 
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(22) 



^f(z,w) ^f(z,w) ^ ^0 ~ a''^'^'^ . . e'y^ - e'^^ - e'^y - (23) 

^xx ^yy a^z.w) ^ ' 

Since deformation potentials and elastic stiffness con- 
stants for WZ InP are not yet available in the literature, 
, , 1 , , we will consider here that the strain effect appears only 

^zz — ^{iii)^xx \^'^) in the ZB structure. This assumption is not totally un- 

realistic because WZ is the dominant phase in the NW. 
Figure [6] shows the effect of strain at /c = for the diag- 
2qM onalized Hamiltonian (k-p and strain terms) as a function 

e'J^) = l^A'^) Q4) of the percentage of strain tensor. A linear variation for 

zz ^ ^ 

^33 the conduction band and the HH band is observed, how- 

ever, the LH and SO bands have a non-linear behavior. 
The order of the HH and LH bands changes when strain 
^z'"^^ = ^'zx""^ = ^'xy""^ = Q (^^) distensive. Table M lists the ZB parameters used in 
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Figure 6. Strain effect in the band edges of ZB InP as a func- 
tion of the percentage of strain tensor. EL and HH shows a 
Unear variation while LH and SO shows a nonhnear behavior. 
The strain effect removes the HH and LH degeneracy. 



the calculations. 



Table H. ZB InP strain parameters. 



Parameter 



ZB InP 



Deformation potentials (eV) 





-6.3735 


D2 


2.2868 


Ds 


8.6603 




-4.3301 




3.8868 







Dq 







-0.4434 




-6.0 


ac± 


-6.0 



Elastic stiffness constant (GPa) 

Cii 1011 

Ci2 561 

Cm 456 



C. Spontaneous and Piezoelectric Polarization 



Piezoelectric polarization appears when a crystal is 
subjected to strain. In ZB semiconductors grown along 
the [111] direction, the magnitude of the piezoelectric po- 
larization, in the unprimed coordinate system of Fig. IH 
is given by |4l| : 



Pi — 2ei4ejk 



(24) 



where ei4 is the piezoelectric constant for ZB materials, 
(i^j^k) are the cartesian coordinates {x^y^z) in a cyclic 
order and ejk are the strain tensor components. 

Applying the coordinate system rotation in the piezo- 
electric polarization vector components in order to de- 
scribe them in the primed coordinate system we obtain: 



Py - 2P,) = 



P'^ = l={P, + Py + P,) = V3P (25) 



The resulting piezoelectric polarization alongside the 
growth direction is then: 



(26) 



r(lll) 



The spontaneous polarization effect in WZ structure 
is due to the relative displacement between the cations 
and anions when the ratio c/a is different from the ideal 
value in the WZ structure. 

In a heterostructure, the effect of the different polariza- 
tions in each region creates an electric field through the 
whole structure. The net electric field in a determined 
layer, z, due to spontaneous and piezoelectric polariza- 
tions in the system is given by (42 



N 



J = l 

^ / 



(27) 



where j sums all over the layers in the heterostructure 
with polarization P, dielectric constant e and length 1. 
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D. Effective mass equation in reciprocal space 

The envelope function approximation [22|,|43j is applied 
to couple the different crystal structures alongside the 
growth direction in the NW. In each region the wave 
function is expanded in terms of the Bloch functions of 
the corresponding poly type. Thus, the wave function of 
the whole system is given by: 



i(k-r) {WZ.ZB)^^^ 



(28) 



where gi{f) are the envelope functions of the l-Xh basis 
state. 

Considering different Bloch functions for each region, 
the Hamiltonian parameters vary alongside the growth 
direction, making it possible to use the common k-p and 
strain matrices, ([5]) and ([9]), for both crystal structures. 
Moreover, since each crystal structure dictates its sym- 
metry to their respective Bloch functions, some matrix 
elements can be forbidden by symmetry in the region of 
a certain crystalline phase. For example, the Az param- 
eter is zero in WZ regions whereas the Ai parameter is 
zero in ZB regions. 

To represent the growth dependence of the Hamilto- 
nian parameters and envelope functions, the plane wave 
expansion is used. This formalism considers the periodic- 
ity of the whole system allowing the expansion of growth 
dependent functions in Fourier coefficients: 



K 



(29) 



where U ^ are the Fourier coefficients of the function U (r) 
and i? is a reciprocal lattice vector. The Fourier expan- 
sion also induces the change k k-\-K in the k-p matrix. 



IV. RESULTS AND DISCUSSION 

The NW system chosen to apply our model is a 
WZ/ZB/WZ single weh structure. Although a real NW 
is composed of multiple polytypical quantum wells with 
different sizes, the analysis of just a single well can bring 
out the physics of the polytypical interface. The effects of 
lateral confinement are neglected in a first approach, as- 
suming NWs with large lateral dimensions. Strain, piezo- 
electric and spontaneous polarization are also included in 
the single well system. 

When both crystal structures are put side by side, a 
band-offset is created at the interface, originating a con- 
finement profile. The band mismatch is also taken from 
reference [28[. Figure El exhibits the WZ/ZB interface 



for InP in two different schemes: in the left, the ener- 
gies in /c = for the diagonalized Hamiltonian and in 
the right, the diagonal terms of the Hamiltonian. Al- 
though the composition of the states in the diagonalized 
energy bands are the same in ^ = 0, the matrix is not 
constructed in this basis. Since the variation along the 
growth direction of the matrix elements is well defined, 
the scheme in the right is more convenient to analyze the 
potential profile of the system. 
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Figure 7. Left side: Band edge energies at /c = for the 
diagonalized Hamiltonian. Right side: Diagonal terms of the 
Hamiltonian at /c = 0. 



In all performed calculations, the entire length of the 
system is set to 500 A and the width of the ZB region, 

is variable. Figure [8] shows the single well potentials 
with and without the effect of strain. Since the potential 
profiles exhibits a type-II behavior, we expect a spatial 
separation of the carriers: electrons are more likely to be 
in the ZB region and the holes in the WZ region. 

The strain considered here, —0.8%, is a intermedi- 
ate value between two data available in the literature 
from ab initio calculations. Reference [34] shows that 
the deviation between the lattice constant of ZB[111] 
and WZ[0001] is -1.3% and reference [M] shows -0.3%. 
Also, reference [32] suggests a difference slighter than 
0.5% between the lattice constants of the two poly types. 
The effect of strain shallows the potential wells in the 
conduction and valence bands, reducing the confinement 
of the carriers. We expect to have less confined states for 
the strained potential compared to the unstrained one. 

The conduction and valence band structures for the 
potential profile without strain are presented in Fig. [9] 
for three different widths of the ZB region. The calcula- 
tions were performed up to 10% in the T — T direction 
and 100% in the T — A direction. For the valence band 64 
energy states are presented while only 18 are presented 
for the conduction band. Since the system has no asym- 
metric potential, the energy bands are two-fold degener- 
ate in spin, therefore 32 states are visible in the valence 
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Figure 8. Diagonal potential profile of the Hamiltonian for 
the polytypical InP system with and without strain. The 
width of the ZB region, /, can change but the whole system's 
dimension remains constant with 500 A. 




Figure 9. Band structure for the unstrained profile in Figure 
m The solid horizontal line is the top energy of conduction 
band well and the dashed horizontal line is the bottom energy 
of the conduction band well. The F — T direction refers to kx 
and r — A to kz. 



band and 9 in the conduction band. For the three dif- 
ferent values of the conduction energy bands are nom- 
inated, from bottom to top, as EL1-EL9, composed of 
|c7,8) states. The valence bands, from top to bottom, are 
nominated as HHl-19, LHl-4, HH20-21, LH5-7, HH22- 
23, LH8-9 for / = 100 A; HHl-16, LHl-2, HH17, LH3, 
HH18, LH4, HH19, LH5-7, HH20-21, LH8-9, HH22-23 
for / = 160A and HHl-14, LHl-2, HH15, LH3, HH16- 
17, LH4-5, HH18, LH6, HH19, LH7, HH20, LH8, HH21, 
LH9, HH22, LHIO for / = 200 A. 

Since the highest valence band states are HH there is 
no significant anticrossing among the energy bands in 
the F — T direction. The anticrossing is characteristic of 
interactions between HH and LH bands in ZB and WZ 
quantum well structures. A slight anticrossing, however, 
can be seen in the energy region just above — 75meV, 
which is next to the interaction region of the |ci^4) and 
|c2,5) profiles. 

Increasing the value of / we find that the number of 
confined states in the conduction band increases. On the 
other hand, for the valence band the number of confined 
states decreases because the WZ region's width also de- 
creases. 

For the strained potential profile, the band structures 



for three different widths of the ZB region are displayed 
in Figure [TOl The calculations were performed consid- 
ering the same extension for the FBZ of the unstrained 
band structure. For the three different values of the 
conduction energy bands are nominated, from bottom to 
top, as EL1-EL9, composed of \c7^s) states. The valence 
bands, from top to bottom, are nominated as HHl-22, 
LHl, HH23, LH2-7, HH24-25 for / = 100 A; HHl-21, 
LHl-3, HH22-23, LH4-8, HH24 for / = 160A and HHl- 
21, LHl-4, HH22-23, LH5-9 for / = 200 A. 

As expected, the number of confined states for the 
strained profile compared to the unstrained one is 
smaller. Nonetheless, the similar confinement trend is 
visible here when the value of / increases: the number 
confined states in the conduction band increases while in 
valence band decreases. 

An interesting feature presented in the strained band 
structure is the presence of some confined states below 
the top region, around —62 meV. This suggests a confine- 
ment in the intermediate region of the |c2,5) and |c3,6) 
profiles. Note that the coupling of these two profiles 
at ^ = happens because of the off-diagonal spin-orbit 
term. 
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The composition of the energy states at ^ = in the 
band structure is similar for the strained and unstrained 
cases: they are just HH or LH states. There is no CH 
states in the energy range considered here. The major 
contribution for CH states comes from the \cs^q) profile, 
which is the lowest one in both cases. 

The information of the energy states' composition can 
reveal important trends in the luminescence spectra for 
this kind of system. For example, at /c = the dom- 
inant symmetry of the energy states belongs to (x,?/), 
which means that the luminescence spectra is more in- 
tense perpendicular to the grow th direction. However, 
experimental measurements [31] indicates that the inten- 
sities perpendicular and parallel to the growth direction 
are almost similar. Therefore, we expect that the con- 
tribution for the parallel luminescence comes from the 
states at /c 7^ 0. For k points away from the F-point, 
there is a stronger mixing between all the basis states. 
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Figure 10. Band structure for the strained profile in Fig. [51 
The solid and dashed lines have the same meaning as in Figure 
[9l The F — T direction refers to kx and T — A to kz- 



The effect of the ZB region width for both strained 
and unstrained potential profile for the conduction and 
valence band states at /c = is presented in Figure 
[TTl It is possible to observe that the number of con- 
fined states in the conduction band increases as the value 
of / increases. On the other hand, the number of con- 
fined states in the valence band decreases. Nevertheless, 



the effect of the variation of / is more significant for the 
conduction band since the electron effective mass in ZB 
{ml /mo = 0.0795) is smaller than the heavy hole mass 
of WZ (m|^^/mo = 1.273 and m^^/mo = 0.158). The 
same trend is also observed in the strained case. Also, 
since the bottom of the well in the conduction band has 
a higher value in the strained case, we can expect the 
interband transition energies to be blueshifted with the 
inclusion of strain effects. 
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Figure 11. The first 5 states of the conduction and valence 
bands at /c = as a function of the ZB region width /. The 
solid line indicates the top of the conduction band well and 
the dashed line indicates the bottom. 



The presence of strain effects gives rise to the piezoelec- 
tric polarization. For the InP ZB, the value for the piezo- 
electric constant used was ei4 = 0.035 C/m^, taken from 
Ref. [44]. For both crystal structures, the value used 
for static dielectric constant was 12.5. The unknown pa- 
rameter is the spontaneous polarization for WZ InP. Ref. 
[32! suggests that this value is smaller than that of InN 
(— 0.03(7/m^). In an attempt to estimate this value for 
InP, we performed the band structure calculations con- 
sidering a range of values for Pgp. 

The energy of the first 5 conduction and valence band 
states at /c = as a function of spontaneous polariza- 
tion in WZ InP for three different ZB region widths 
is presented in Figure [121 The considered values for 
spontaneous polarization are —0.02 C/m^, —0.015 C/m^, 
-0.01 C/m^, -0.005 C/m^ and -0.001 C/m^. For / = 
160 A and / = 200 A there is a crossing between the con- 
duction and valence band states. This is not observed 
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experimentally therefore we consider this region forbid- 
den. Then, the allowed values for spontaneous polariza- 
tion considered here are then —0.01 C /m?^ —0.005 C /m? 
and -0.001 C/m^. 
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Figure 12. The first 5 states of the conduction and valence 
bands at /c = as a function of WZ spontaneous polarization 
Psp. Notice the crossing of valence and conduction bands. 
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Figure 13. Diagonal potential profile of the Hamiltonian for 
the polytypical InP system considering strain, piezoelectric 
polarization in the ZB region and spontaneous polarization in 
the WZ region for / = 200 A. The values for the spontaneous 
polarization were chosen from Figure [T2l 



The diagonal potential profile including effects of 
piezoelectric and spontaneous polarization is presented 
in Figure [131 The ZB region is fixed at 200 A. Analyzing 
these profiles, we expect to have strong coupling in the 
band structure for higher values of Pgp since the profiles 
are more close to each other. This induces the mixing of 
states because an energy value can be in more than one 
profile. 

The resulting band structures for the three differ- 
ent potential profiles of Figure [13] are shown in Figure 
[T4l For the three different values of Pgp, the conduc- 
tion energy bands are nominated, from bottom to top, 
as EL1-EL9, composed of \c7^s) states. The valence 
bands, from top to bottom, are nominated as HHl-2, 
LHl, HH3, LH2, HH4-5, LH3, HH6, LH4, HH7, LH5, 
HH8, LH6, HH9, LH7, HHIO, LH8, HHll, LH9, HH12, 
LHIO, HH13, LHll, HHM, LH12, HHIS, LH13, HH16, 



LH14, HH17, LH15 for Psp = -0.01 C/m^; HHl-3, LHl, 
HH4, LH2, HH5-6, LH3, HH7, LH4, HH8, LH5, HH9- 
10, LH6, HHll, LH7, HH12, LH8, HH13, LH9, HHM, 
LHIO, HH15-16, LHll, HH17, LH12-13, HH18, LH14 
for Psp = -0.005 C/m^ and HHl-7, LHl, HH8-9, LH2, 
HHlO-11, LH3, HH12-13, LH4, HH14-15, LH5, HH16, 
LH6, HH17-18, LH7, HH19, LH8, HH20, LH9, HH21-22, 
LHIO for Psp = -0.001 C/m^. The number of HH states 
increases as the value of Psp decreases. 

The anticrossings and also the spin splitting in the va- 
lence sub bands are more visible for higher values of Psp. 
The strength of the resulting electric field not only in- 
creases the mixing of HH and LH states but also increases 
the value of the spin splitting in each sub band. This spin 
splitting is known as the Rashba effect [45] and is due 
to potential inversion asymmetry, even though the term 
a{a X k) • E does not appear explicitly in the Hamiltonian 



12 



The number of confined states decreases as the sponta- 
neous polarization decreases. On the other hand, the en- 
ergy difference between the conduction and valence band 
ground state increases as the spontaneous polarization 
decreases, blueshifting the interband energy transitions. 




Figure 14. Band structure for the profiles presented in Figure 
1131 The spin-splitting of the energy bands is due to the field 
induced asymmetry. 



The effect of piezoelectric and spontaneous polariza- 
tion also induce carriers' spatial separation. This effect, 
in the probability densities in A? = can be seen in Fig. 
[T5l The lowest four states of the conduction band and 
the highest four states of the valence band are presented. 
At /c = the wave functions of spin-up and spin-down 
are degenerated. We can see that the overlap increases 
for more excited states, also blueshifting the energy peak 
in the interband transitions. Since the potential profile 
is not completely even or odd, the envelope functions no 
longer have well defined parities. 



V. CONCLUSIONS 

The basic result of this study is the theoretical model 
based on the k-p method and group theory concepts to 
calculate band structures of WZ/ZB polytypical systems 
in the vicinity of the band edge. The method allows us to 
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Figure 15. Probability densities at /c = for the lowest four 
states of the conduction band and the highest four states of 
the valence band of Fig. 1141 The solid lines are the HH states 
and dashed lines represent the LH states. 



describe in the same matrix Hamiltonian the ZB and WZ 
structures, with kz along the [111] and [0001] directions, 
respectively. Since the WZ structure is less symmetric, 
the ZB parameters are assigned to the WZ ones. Our 
method not only is able to describe the k-p terms of the 
Hamiltonian but also includes the strain and polarization 
(spontaneous and piezoelectric) effects. 

Extracting the parameters of WZ InP from Ref. jlH 
we applied our model to a WZ/ZB/WZ single well in 
order to understand the physics of the polytypical inter- 
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face. The potential profile at the interface WZ/ZB is 
type-II, whose feature is the spatial separation of carri- 
ers. The performed calculations in this study holds this 
characteristic. 

Due to the lack of parameters in the literature for WZ 
InP, only the strain effect in the ZB region was consid- 
ered here. This seems to be a reasonable consideration 
since the WZ structure is the dominant phase in NWs 
structures. However, such strain parameters would be 
fundamental in a system that the stable lattice constant 
is a intermediate value between WZ and ZB InP lattice 
parameters. 

Within the limitation of strain, the piezoelectric polar- 
ization was also considered in the ZB region. For the WZ 
region, only the spontaneous polarization appears. Since 
there is no value in the literature for the spontaneous po- 
larization of WZ InP, a range of values were considered in 
the simulations. Some of these values, however, induces 
a negative gap in the system. There is no data in the 
literature that corroborates this effect. 

The proposed model, jointly with the obtained results, 
proved to be useful in the study of electronic band struc- 
tures of WZ/ZB poly typical systems, such as NWs. Ex- 
ploring the opportunities of band gap engineering con- 
sidering not only different compounds, but also different 
crystal structures, could lead to the development of novel 
nano devices. 
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